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Abstract. Since a large fraction of active galactic nuclei (AGN) is missed in common UV-excess surveys and is even hard to 
find in radio, near-IR and X-ray surveys, we have used a new AGN selection technique which is expected to be not affected 
by extinction. Within the scientific verification of the ISOCAM Parallel Survey at 6.7 /im we have discovered objects with 
exceptional mid-infrared (MIR) emission. They are essentially not detected on IRAS-ADDSCANs and only very few of them 
show up in the NVSS and FIRST radio surveys. Various colour criteria of the 6.7 fim data with 2MASS and optical wavebands 
show that the sources reach more extreme IR colours than the sources in the Hubble Deep Field-South and the ELAIS survey. 
The comparison with known object types suggests that we have found AGN with a pronounced MIR emission, probably due to 
circum-nuclear dust. First results from optical spectroscopy of ten candidates corroborate this interpretation showing four AGN, 
two reddened LINER and four extremely reddened emission-line galaxies with MIR/FIR flux ratios higher than for known pure 
starburst galaxies. The results will make a significant contribution to the debate on the entire AGN population. 
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• • i 1. Introduction pie and concluded that up to 80% of the quasars have been 

missed in conventional optical surveys, provided that the red- 



The realisation that obscuration plays a critical role in the clas- der colours of the radio _i oud quasars are due t0 dust redd en 
^ sification of AGN fundamentally inspired the current research. ing and that me radio . quiet quasars contain as much dust as the 
Ci ; Attempts to overcome the limitations of dust extinction and radio . loud ones . By new strateg ies in the optical, assuming that 
to identify the entire AGN population - including type 2 and the narrow . Une regions ^ sufficiently extended and that only 
dust-enshrouded AGN - encompass surveys in the near-IR, ra- the continuum emission is hidden, type-2 quasar candidates 
dio, and X-ray regimes. However, searching for very red AGN have been selected afj objectfj wkh nam)W permitted emission 
the colour selection via J - K B > 2 (Cutn et al. 2001) ex- lines an d high [OIII]A5007 equivalent widths (Djorgovski et 
eludes most of the known AGN (Barkhouse and Hall 2001), al 20 oi, Zakamska et al. 2003). Applying a moderate colour 
only about 30% of AGN are radio-loud (Urry & Padovani cut j _ Kg > l % from about 1500 sources in 2MASS Francis 
1995), and there seems to exist many X-ray faint AGN (Wilkes et al (2004) find only tentat i ve evidence that Seyfert2 nuclei 
et al. 2002). Thus, a considerable fraction of the AGN popu- m m common in the NIR selected survey than in blue se- 
lation must have escaped detection due to observational bias. lected galaxy surveys , and the y can place only very weak con- 
Webster et al. (1995) found that their radio-selected quasar straints on any popu i ation f dusty AGN. 
sample is significantly redder than an optical comparison sam- 

~Z , rr 7Z ~ 7". n n. <a . u a \ The disadvantage of heavy extinction in optical and NIR 

Sena offprint requests to: Martin Haas (haas@astro.rub.de) ° J * 

* Based on observations with the Infrared Space Observatory ISO, surve y s can turn into a valuable detection tool, when observing 

an ESA project with instruments funded by ESA Member States (es- dust-surrounded AGN at MIR wavelengths. There, the reemis- 

pecially the PI countries: France, Germany, the Netherlands and UK) sion of the hiding dust heated by the strong radiation field of the 

and with the participation of ISAS and NASA, and also based on ob- AGN should be seen easily as MIR excess. In fact, for known 

servations at the South African Astronomical Observatory SAAO. (powerful) AGN of both type 1 and type 2 a steep near- to mid- 
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IR slope has been revealed by sensitive MIR observations (e.g. 
Haas et al. 2003, Haas et al. 2004, Siebenmorgen et al. 2004). 

We therefore started a new approach, searching for AGN 
by means of their MIR emission of the nuclear dust torus. 
However, one complication with this method has to be solved: 
Since luminous IR starburst galaxies may also show a pro- 
nounced MIR emission due to the PAH emission bands around 
7.7 /im, it is of special importance to distinguish them from 
AGN. In this Letter we describe the new technique for the mid- 
IR selection of AGN candidates using IR colour diagrams and 
report about first results from optical spectroscopy. 

2. Selection of MIR sources 

ISO has performed a serendipitous survey at 6.7 /im (LW2 
band), the ISOCAM Parallel Survey, with 6" spatial resolution 
and a positional accuracy of better than 3" (Siebenmorgen et 
al. 1996, Ott et al. 2003). Over 27 square degrees of the sky are 
processed and currently being catalogued. For point sources the 
3 a detection limit is about 0.5 mJy. 

Within the scientific verification of the 17000 detected 
sources (Ott et al. 2004 in prep.) we have selected unresolved 
sources at galactic latitude |6| > 20°. We then performed cross 
correlations with the 2MASS all sky point source catalogue 
(Cutri et al. 2003), with the USNO-B, DSS and UCAC op- 
tical catalogues, as well as the NVSS and FIRST radio sur- 
veys, and analysed IRAS ADDSCANs. We excluded objects 
which have multiple NIR and optical counterparts within 10", 
or are contaminated by extended sources (2MASS XSC), or 
have proper motion (pm > 3 a from UCAC). The resulting list 
contains about 3000 ISOCAM point sources with BRJHK 
and LW2 (6.7 /im) photometry (henceforth denoted ISOCP 
sources). The brightness ranges are: B — 14.5 . . . 20.5 mag, 
K = 12.5 ... 16 mag, LW2 = 8 ... 12 mag ±0.2 ... 0.4 (Vega 
based system). The typical errors on the H—K s and K S — LW2 
colours are 0.1 ... 0.2 mag and 0.3 . . . 0.5 mag, respectively. 

3. AGN candidate selection 

The ISOCP list contains various object classes like stars, nor- 
mal and active galaxies. Fig. [2 (top) shows their distribution 
in the H — K s versus K s — LW2 colour-colour diagram. 
Noteworthy, the ISOCP sources clearly reach more extreme IR 
colours than those of the Hubble Deep Field South (Oliver et al. 
2002, Mann et al. 2002) and the ELAIS fields (Vaisanen et al. 
2002, Rowan-Robinson et al. 2003). The corresponding limits 
are indicated in Fig.^ 

To identify the types of objects that populate different re- 
gions of the colour-colour diagram shown in Fig. 1 (top), we 
show the same diagram for sources with known identification 
in Fig. 1 (bottom). Stars, i.e. proper motion objects, lie in the 
lower left corner at K s - LW2 < 1 and H - K s < 0.5. This 
area also contains radio-quiet elliptical galaxies from the re- 
vised Shapley-Ames catalog. 

UV-excess selected Palomar-Green (PG) quasars (AGN 
type 1, data from Haas et al. 2003) and powerful FR2 radio 
galaxies from the 3CR catalogue (AGN type 2, data from Haas 
et al. 2004 and Siebenmorgen et al. 2004) lie in the upper right 
corner at K s — LW2 > 2.7 and H — K s > 0.5 - a region 
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Fig. 1. Top: H — K s vs. K s — LW2 colour-colour diagram 
of the 3000 unresolved ISOCAM sources at |fe| > 20°. Circles 
and fat dots denote the AGN candidates, small dots all other 
sources. The dotted oblique lines indicate the limits reached by 
the HDF-S and ELAIS sources, lying below and left hand of 
the lines. Bottom: Comparison samples. 

that is also populated by 3CR quasars (AGN type 1). As an ad- 
ditional check, we cross-correlated the UV-selected quasars in 
the Veron-Cetty & Veron catalog (2003, 1 1th edition) with the 
2MASS catalog. More than 90% of the quasars with V < 18 
mag and z < 0.8 lie at H — K B > 0.5. At z ~ 1.5 the 
fraction goes down to 10%, probably due to the shift of the 
NIR excess beyond the iT s -band. The sources from the PG and 
the 3CR catalogs, respectively, shown in Fig.[Jbottom are at 
z < 0.8. For simplicity, we denote the sources in the colour 
range K s — LW2 > 2.7 and H — K s > 0.5 as excess sources. 

The ISOCP excess sources are neither carbon stars (e.g. 
Liebert et al. 2000) nor AGB stars (Joris Blommaert, priv. 
comm.), both types showing optical and NIR colours different 
from these ISOCP sources. Also, dusty young stars show dif- 
ferent optical colours and they are not expected at high galactic 
latitude |6| > 20°. 

Normal star forming galaxies (Boselli et al. 1998, Dale 
et al. 2001) populate the full range of K s — LW2, reach- 
ing into the regime of AGN, but they lie significantly lower 
in H — K s (< 0.5), so that they can be distinguished from 
the PG/3CR AGN. Even IR luminous starburst galaxies hav- 
ing 10 11 < Lpir [L ] < 10 12 , like Arp 157, NGC 1569, 
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Table 1. Mean colours of Seyfert galaxies from Clavel et al. 
(2001) and ULIRGs from Klaas et al. (2001). N gives the num- 
ber of objects used for deriving the two colours, respectively. 
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NGC2146, NGC3256 or M82, fall into this colour range 
(H -K s < 0.5). 

Seyfert 1 and 2 galaxies, from the sample by Clavel et 
al. (2001) with averages listed in Tab.[Q lie in the "transi- 
tion range" between normal star forming galaxies and PG/3CR 
AGN. While Sy 2s on average are below H — K s = 0.5, Sy Is 
are above and for clarity of the plot only the word "Seyferts" 
delineates their colour location in Fig.[2bottom. 

The colours of IR ultra-luminous (Lfir > 10 12 L ) star- 
burst galaxies (including LINERs), which we derived from 
the sample of Klaas et al. (2001) together with 2MASS and 
ISOCAM 6.7/itn data, he in the range K s - LW2 > 2.7 and 

0. 2 < H - K s < 0.8 (Tab.GJ. 13 of them lie below and 8 above 
H-K s = 0.5 overlapping with the range of the PG/3CR AGN. 
In Fig.Honly the word "SB-ULIRGs" is placed at their mean 
location. On the other hand, the AGN-ULIRGs with "warm" 
^12/^25 colours (like Mrk231 and IRAS 05189-2524) clearly 
fall into the colour range of the PG/3CR AGN (Tab.[T}. 

In following we discuss, whether the ISOCP excess sources 
could be explained in terms of pure starburst IR galaxies. We 
will argue that such starforming galaxies would be either too 
bright or too distant to match both the observed fluxes at K s , 
LW2 and at 60/im and the corresponding colours of the ISOCP 
excess sources. Our argumentation is essentially based upon 
the following empirical facts: 

1. The ISOCP excess sources have median (mean) fluxes 
at K s and LW2 of 0.66 mJy (0.9 mJy) and 1.75 mJy 
(3.2 mJy), respectively. 

2. The LW2 flux of a starburst galaxy with z < 0.15 is mostly 
due to the PAH bands around 7.7 /im (e.g. Laurent et al. 
2000). For z > 0.15 this feature is shifted out of the LW2 
passband (5 — 8.5 /im) resulting in a reduced K s — LW2 
colour below 2.7. 

3. If a red (H — K s > 0.5) starburst galaxy has a K s — band 
excess - originating from hot dust close to the sublima- 
tion temperature - then its 60 /im emission is also ex- 
pected to be strong. E.g. the five SB-ULIRGs observed by 
Klaas et al. (2001) with H - K s > 0.4 show a flux ra- 
tio Fqq/F^^ ~ 230 . . . 570. Less luminous, more gently 
star-forming galaxies might have a smaller Fqq/ Flw2 ra- 
tio, but - on the other hand - the typical known SB-LIRGs 
like M 82 are bluer (H-K s < 0.5) than the present ISOCP 
excess sources. 

PAH bands: First, we make use of the PAH bands and dis- 
cuss the ISOCP excess sources with F2.2fj.rn 1 m Jy- In order 
to match this limit, a typical SB-ULIRG, like e.g. IRAS 17208- 
0014, with F 2 .2 M m « 20mJy at z = 0.042, would have to 
lie at z = 0.188. In this case, however, K s — LW2 would 



attain values below 2.7, see point (2.) above. The same argu- 
ment applies to other SB-ULIRGs (like MRK273, ARP220, 
IRAS 14348-1447 and IRAS 23365+3604). Hence, a pure star- 
burst counterpart of the faint ISOCP excess sources must have 
a lower luminosity, i.e. being at most an SB-LIRG. 

-^60 / Flw2 ratio: Second, we make use of the F 60 /Flw2 
ratio, which attains typically values above 230 for SB-ULIRGs. 
A similar flux ratio is expected for highly dust-enshrouded SB- 
LIRGs with H - K s > 0.5 (down-sized SB-ULIRGs), if they 
exist. In order to match Flw2 > 1 mJy for the bright ISOCP 
sources, the expected flux at 60 /im would have to be at least 
F 60 230. . .570mJy, which is above the IRAS detection 
threshold at low-cirrus high-galactic latitudes. 

For faint ISOCP sources with fluxes Flw2 < 1 mJy we 
expect useful IRAS upper limits. We have examined the 70 
unclassified ISOCP excess sources (Fig.Q top) individually: 
only eight of them show a marginal detection in the 60 /im 
IRAS ADDSCAN data, the remaining sources have 3cr up- 
per limits between 90 and 190 mJy. Two of the detected 
sources exhibit Fqq/ Flw2 < 230, all of the remaining sources 
have Fqq/Flw2 < 180 and at least 40 of them even have 
FQ0/FLW2 < 100. These low limits show that most of the 70 
ISOCP sources are characterized by a high MIR/NIR flux ra- 
tio, which is not accompanied by a high FIR/MIR ratio as ex- 
pected for known dust-enshrouded pure starburst galaxies with 
H - K s > 0.5. 

If the ISOCP excess sources are not a new population of 
star-forming galaxies, which are highly dust-enshrouded (with 
red H — K), rich in PAH emission (with bright F^y^) and have 
cool FIR colours (with low 60/xm flux), then the arguments 
above favour a significant AGN contribution in these sources. 
Furthermore, the four AGN-ULIRGs show an Fqq/Flw2 ra- 
tio between 10 and 50, consistent with the 60/im upper limits 
found for the ISOCP excess sources. 

Based on the comparison with known object types we pre- 
dict that most of the 80 ISOCP excess sources house an AGN 
and have Seyfert luminosities, but could be more powerful if at 
z ;> 0.2. Therefore, we consider the ISOCP excess sources as 
promising AGN candidates. At redshift about z > 0.8, when 
the AGN-typical SED bump shifts, the colours of the sources 
may become bluer than H — K s = 0.5, and a refined analysis 
using other filters has to be applied in order to identify all AGN. 
Remarkably, the new method is - a priori - not biased against 
optical-UV selected QSOs. For comparison, the 2MASS AGN 
search catches only the range H — K s > 0.7 (roughly cor- 
responding to the J — K s > 1.2 criterion used by Francis et 
al. 2004), hence may ignore more of the known AGN - even 
among the local ones. 

4. Optical spectroscopy 

Obviously, our predictions about the nature of the ISOCP ex- 
cess sources have to be verified by optical spectroscopy. Ten 
of our 80 AGN candidates are already listed in the NED as 
QSOs. None is listed as a star (in SIMBAD) or starburst galaxy. 
Twelve sources have 1.4 GHz radio detections, 4 being radio 
loud. Only one source (3C345, Sandage & Wyndham 1965) 
has also been identified spectroscopically as QSO in the Sloan 
Survey (SDSS DR1, Abazajian et al. 2003) while SDSS DR2 
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Fig. 2. Optical spectra of three ISOCP sources with H - K s > 0.5 and K s - LW2 > 2.7 obtained at the SAAO 1.9 m 
telescope: ISOCP02831: type-1 AGN with broad H Q and Hp lines, ISOCP08080: type-2 AGN, ISOCP10881: reddened emission 
line galaxy with H a /Hp > 10. 



contains only one more spectrum of an ISOCP excess source 
(RA,Decj 20 oo = 13:25:07.8, +05:41:07)- it is a QSO at z = 0.2. 

Because of the marginal spectroscopic classification of the 
ISOCP sources, we have started a spectroscopic survey for the 
70 unknown ISOCP excess sources (K s — LW2 > 2.7 and 
H - K s > 0.5). First results obtained at the 1.9 m SAAO 
telescope by end of January 2004 on ten sources confirm their 
extragalactic nature (Leipski et al. 2004 in prep.). The obser- 
vations reveal two type-1 AGN and two type-2 AGN, one of 
which (ISOCP08080, Fig.fJ has a high [O III]A5007 luminos- 
ity above 10 8 L and can be classified as QSO-2. Furthermore, 
we find two reddened LINER and four extremely reddened 
emission-line galaxies. Example spectra are given in Fig. [2] 
The visual extinction of the emission-line galaxies as infeiTed 
from the ratio H a fH/3 > 10 is Ay > 3 mag; actually Ay is 
expected to be much stronger when derived from data at longer 
wavelengths, as found in many dust-enshrouded IR sources 
(e.g. Haas et al. 2001). The redshift range is 0.1 < z < 0.5. 
Evidence is growing that LINER galaxies contain an AGN (e.g. 
Satyapal et al. 2004). Also the four emission-line galaxies may 
contain an AGN, since their flux ratio Fqo / Flw2 < 100 is be- 
low that found for known pure starburst IR galaxies. Therefore, 
the ongoing results with a fraction of 40% (4/10) optically iden- 
tifyable classical AGN are consistent with our prediction that 
most of the ISOCP excess sources are AGN. Alternatively, we 
would have to postulate a new population of dust-enshrouded 
unusually cool star-forming galaxies. 



5. Conclusions 

We have discovered a sample of unique MIR excess sources. 
Various arguments suggest that they likely contain an AGN. 
The ongoing optical spectroscopy indicates 40% of the sam- 
ple to be classical AGN and the remaining part to be dust- 
enshrouded sources. They appear enigmatic in terms of having 
a high MIR/NIR flux ratio, but relatively low FIR upper lim- 
its. Since a moderate FIR/MIR flux ratio is more typical for 
AGN, it argues against pure starbursts. This conjecture has to 
be confirmed. Because of the high extinction, optical obser- 
vations might not be able to see the true AGN or starburst or 
composite nature of those sources. Therefore, future observa- 
tions with XMM-Newton, the VLA and IR spectroscopy with 
the Spitzer Space Telescope are planned. 



If all 80 ISOCP MIR excess sources turn out to be AGN 
and if 32 of them (= 40%) are classical AGN which could have 
been identified also in optical surveys, then the number counts 
of AGN may increase by a factor of up to 80/32 = 2.5 due to 
the inclusion of dusty AGN. Clearly, a reliable comparison will 
have to wait until spectra of more ISOCP sources are obtained. 

The emphasis of this letter is on presenting a promising 
strategy how to find the presumed missing AGN. The selec- 
tion of MIR excess sources via the H — K s and K s — LW2 
colour-colour diagrams in combination with the low 60/im up- 
per flux limits turns out to be efficient in separating starbursts 
from AGN. It will establish a long sought-after technique for 
selecting the full population of AGN largely free of extinction. 
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